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ABSTRACT 

A weakly charged flexible polyelectrolyte chain in a neutral spherical cavity is analyzed using 
self-consistent field theory (SCFT) within an explicit solvent model. Assuming the radial 
symmetry for the system, it is found that the confinement of the chain leads to creation of a 
charge density wave along with the development of a potential difference across the center of 
cavity and the surface. We show that the solvent entropy plays an important role in the free 
energy of the confined system. For a given radius of the spherical cavity and fixed charge 
density along the backbone of the chain, solvent and small ion entropies dominate over all 
other contributions when chain lengths are small. However, with the increase in chain length, 
chain conformational entropy and polymer- solvent interaction energy also become important. 
Our calculations reveal that energy due to electrostatic interactions plays a minor role in the 
free energy. Furthermore, we show that the total free energy under spherical confinement 
is not extensive in the number of monomers. Results for the osmotic pressure and mean 
activity coefficient for monovalent salt are presented. We demonstrate that fluctuations 
at one-loop level lower the free energy and corrections to the osmotic pressure and mean 
activity coefficient of the salt are discussed. Finite size corrections are shown to widen the 
range of validity of the fluctuation analysis. 
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I. INTRODUCTION 



Polyelectrolytes are ubiquitous in nature and exhibit rich behavior. In the past, a great 
deal of theoretical^'^'^'^'^'^'^'^'^'^ and experimental effort3^>^>^>^>^»^ have been made in 
understanding their characteristics in solutions with concentrations ranging from dilute to 
concentrated. One of the remarkable discoveries in the last century is the theory of simple 
electrolytes made by Debye and HiickeU^, where electrostatic interactions get screened (col- 
loquially referred to as the Debye screening). Similar phenomenon was shown to be present 
in the case of neutral polymers, where monomers interact by short range excluded volume 
interaction potential (known as Edward's screening^^) . In the case of polyelectrolytes, both 
kinds of screening effects are present and their coupling via the concept of double screening 
was introduced by Muthukumar^, and the behavior of polyelectrolyte solutions was de- 
scribed in terms of these screening phenomena. As a result, different concentration regimes 
in polyelectrolyte solutions were predicted and verified experimentallyii»i^. 

However, most of the computer and real-world experiments involve finite vol- 
ume, where boundary effects play a significant role and can not be ignored. Re- 
cently, there has been a resurgence of interest in studying polymers within confined 
domainsi^i22i2ii2^i2^iMi2^i2^i2Ii22ii22ii^. Unlike neutral polymers, little is known about the physics 
of polyelectrolytes under confinement, a situation realized in liposome- mediated delivery of 
macromolecules to the cells^^, translocation experiments involving RNA/DNA— i^i, synthetic 
polyelectrolytes^ etc. Underlying physics in these experiments is governed by the confine- 
ment effects on a single polyelectrolyte chain. Physically, confinement forces interaction 
among the monomers and the conformational entropy of the chain gets lowered due to less 
number of conformational states available to the chain. If small components like salt ions 
and solvent molecules are also present in addition to the polymer, then translational entropy 
of these components gets diminished, making confinement a thermodynamically unfavorable 
process. Quantitative estimates of these confinement effects are desirable. Recently, a quan- 
titative description of the finite size corrections^ to free energy for electrolytic systems has 
been presented. The analog of these calculations for polyelectrolyte systems has not been 
attempted yet and is one of the goals of this study. 

In this study, we focus on a single polyelectrolyte chain confined in a neutral spherical 
cavity and use radial symmetry to obtain the mean field results. This theoretical model 



is pertinent in understanding many important physical processes. Few promising applica- 
tions of this model are the computation of free energy barriers for the chain to move out of 
the confining space, osmotic pressure of polyelectrolytes, etc. To start with, we consider a 
situation where inner and outer dielectric constants of medium are different (say ei and 
respectively). For this situation, the electrostatic potential is to be calculated using continu- 
ity of electrostatic potential and normal component of displacement vector at the boundary. 
It turns out that these boundary conditions are equivalent to a continuous dielectric medium 
(of dielectric constant everywhere) with charges confined within the boundary of sphere 
for the radially symmetric globally electroneutral system. In other words, dielectric mismatch 
effects disappear due to the use of radial symmetry. We must point out that radial symmetry 
is strictly valid only at the mean field level and fiuctuations break this symmetry. In the 
concentrated regime, we have been able to compute the fiuctuation contributions without 
recourse to radial symmetry. 

Unlike cylindrical^^ and rectangular— i2i confinements, a single self-avoiding chain in a 
spherical cavity has been shown to be a polymer solution problen*^'^'^'^'^'^ with dif- 
ferent degrees of confinements corresponding to different concentration regimes as seen in 
polymer solution theories. We are here interested in estimating different thermodynamic 
contributions to the free energy of confined polyelectrolytic system. Computing these contri- 
butions using simulations is a formidable task. However, self-consistent field theory (SCFT) 
presents a faster and an accurate way to address this problem. SCFT allows us to explore 
the role of confinement in free energy at the mean field level (also known as the saddle 
point approximation) and then capture the role of fiuctuations by expanding free energy 
functional around the mean field solution (one- loop calculations). Unlike earlier studies on 
single neutral chairt^^, we use "explicit solvent model", which captures solvent entropy in 
a more realistic way. At the mean field level, free energy for a self-avoiding chain under 
strong confinement (in vacuum) has been shown to be proportional to wppQ, where w is a 
measure of the strength of excluded volume interactions, pp is monomer density and Q is the 
volume of the cavity. Remarkably, same proportionality is exhibited by polymer solution 
theories^i^ in the concentrated regime. However, at low polymer concentrations (i.e. dilute 
and semi-dilute regime), fiuctuations around the saddle point solution become important 
and saddle point approximation breaks down. In that case, non-perturbative techniques 
have to be devised to compute the correct free energy. In this work, we explore the free 



energy in the concentrated regime, when radius of gyration of the chain is comparable to 
the radius of the cavity and mean field theory is still applicable. 

Dividing the free energy into energetic and entropic parts using thermodynamic argu- 
ments, we have identified the role of individual components. Moreover, one loop calcula- 
tions provide insight about the corrections to the bulk expressions from fluctuation effects 
in the concentrated regime. For a globally neutral system with only small ions (i.e. without 
polyelectrolyte) inside a neutral cavity, local electroneutrality is the equilibrium state. In 
contrast, local electroneutrality is violated in the presence of a polyelectrolyte chain due to 
the depletion effects present in the system. We have studied the resulting monomer and 
charge density distribution for different sets of relevant parameters of the problem. To link 
with the experiments, we have computed the osmotic pressure and mean activity coefficient 
for monovalent salt. 

This paper is organized as follows: theory is presented in Sec. [Tll calculated results and 
conclusions are presented in Sec. [ITT] and WV\ respectively. 

II. THEORY 

A. Self-consistent field theory 

We consider a spherical cavity of radius R containing a single flexible polyelectrolyte chain 
of total Kuhn segments, each with length b. The polyelectrolyte chain is represented as 
a continuous curve of length Nb and an arc length variable t is used to represent any seg- 
ment along the backbone. We assume that there are monovalent counterions (positively 
charged) released by the chain (and assuming that the chain is negatively charged for the 
sake of specificity). In addition, there are ions of species 7 (= +, — ) coming from added 
salt (in volume Q = AtcR^/S ) so that the whole system is globally electroneutral. Let Zj be 
the valency (with sign) of the charged species of type j. Moreover, we assume that there are 
Hg solvent molecules (satisfying the incompressibility constraint) present in the cavity and 
for simplicity, each solvent molecule occupies a volume (vg) same as that of the monomer 
(i.e. Vg = b^)- Subscripts p,s,c,+ and — are used to represent monomer, solvent , counte- 
rion from polyelectrolyte, positive and negative salt ions, respectively. Degree of ionization 
of the chain is taken to be a and we consider smeared charge distribution so that each of 



the segments carries a charge eaZp^ where e is the electronic charge. 

If electrostatic interactions between solvent molecules and small ions are ignored, the 
partition function for the system can be written as 
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where R(t) represents the position vector for t*^ segment and subscripts j, a = s, c, +,— . 
In the above expression, it is understood that the factor of 1/2 in the last term inside 
the exponent is there, only when j = a. In Eq. ( [1]), ksT is the Boltzmann constant 
times absolute temperature. Vpp{r),Vss{r) andV^s(r) represent monomer-monomer, solvent- 
solvent and monomer-solvent interaction energies, respectively, when the interacting species 
are separated by distance r =| r | and are given by 

Vpp{r) = WppSiv) + (2) 
Ks(r) = WssS{r) (3) 
Vpsir) = Wps5{r). (4) 

Here, Wpp,Wss and the excluded volume parameters, which characterize the 

short range excluded volume interactions between monomer-monomer, solvent-solvent and 
monomer-solvent, respectively. 5(r) is the three dimensional Dirac delta function and e is 
position independent effective dielectric constant of the medium ( in units of 47reo, where 
€0 is the permittivity of vacuum). Delta functions involving microscopic densities enforce 
incompressibility condition at all points in the system (po being total number density of the 
system). Interactions between polyelectrolyte segments and small ions as represented by 
Vpj, are given by 

Vp,{v) = W - for j = c,+,-. 5 

In writing the interaction energies between polyelectrolyte segments and small ions, we 
have taken small ions to be point charges so that they have zero excluded volume and 



interactions are purely electrostatic in nature. In the point charge limit for the small ions, 
interaction between two small ions is written as 



^^■"^""^ " ^fc^T r ^^'^ J,a = c,+,-. (6) 
Following the standard protocoF^ to obtain saddle point equations (see Appendix A), the 
Poisson-Boltzmann equation for electric potential gets coupled to the well-known modified 
diffusion equation for chain connectivity. In particular, the saddle point equations are given 

by 
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These equations are equivalent to those derived by Shi and Noolandi^, and Wang^^ et 
ai, although the method of derivation is different as briefly outlined in Appendix A. In 
these equations, P/3(r) and 0/3 (r) are respectively the macroscopic number density and the 
field experienced by particles of type P, due to excluded volume interactions. All charged 
species experience electrostatic potential represented by tp{r) above. Note that 4j{r) in 
above equations is dimensionless (in units of ksT/e) and Ib depicts the Bjerrum length 
defined as = e^/ekBT. Moreover, Xps is the dimensionless Flory's chi parameter and 
Po = {N + Us)/^ = Also, 7]{r) is the well-known pressure field introduced to enforce 

the incompressibility constraint. The function g(r, s) is the probability of finding segment 
s at location r, when starting end of the chain can be anywhere in space, satisfying the 
modified diffusion equation^ 
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Also, Qp represents the partition function for the particle of type j3 in the field experienced 
by it, given by 
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Using the above equations, approximated free energy at the extremum (saddle point ap 
proximation) is given by 
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where j = s, c, +, — and F^/kBT = y {^Wpp + UsWss) is the self-energy contribution aris- 
ing from excluded volume interactions. Using thermodynamic arguments^ and assuming 
dielectric constant (e) to be independent of temperature (T), the free energy (Eq. ( fT8l) ) is 
divided into enthalpic contributions due to excluded volume, electrostatic interactions and 
entropic part due to small ions, solvent molecules and the polyelectrolyte chain. Denoting 
these contributions by E^, E^, Sions, Sgoivent and Spoiy, respectively, the free energy is written 
as 
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So far, we have presented a general field theoretical treatment for a single polyelectrolyte 
chain and haven't considered confinement. We study the role of confinement by solving Eqs. 
( [3- [T^ under a particular set of boundary conditions and constraints, which are presented 
in the next section. Also, the limits of volume integral in Eqs. ([19]- [211) vary over the volume 
of confining spherical cavity. 



B. Boundary Conditions and Constraints 

The above treatment leads to coupling of non-linear Poisson-Boltzmann equation with 
modified diffusion equation. Both of these equations are second order differential equations 
and hence, two conditions are required for each, in order to obtain a unique solution. In 
addition, for s dependent diffusion equation, an initial condition is needed to start the 
computations. To solve these equations, we exploit the assumed spherical symmetry of the 
system so that g(r, s) q{r, s) and the Laplacian is given by 

V; = ~ {r'§;) (25) 

Due to symmetry of the system, additional requirements need to be fulfilled by the solution. 
Here, we summarize all these conditions: 



Boundary Conditions : — - — \r=R= 0, q(R, s) = for all s (26) 

or 

Initial Conditions : g(r, 0) = 1 for all r 7^ i? (27) 

dipiv^ Oq(t^ 5) 

Symmetry Conditions: — — |r=o= — tt — |r=o= for all s (28) 

or Or 

Boundary and initial conditions for q{r, s) correspond to the facts that (a) segments are 
excluded from boundary so probability of finding any segment at the boundary is zero and 
(b) the ends can be anywhere inside sphere. Symmetry condition for q{r, s) is invoked 
because we are looking for a symmetrical solution of monomer density about the center 
without any discontinuity. Boundary condition and symmetry conditions for electrostatic 
potential are obtained by using the Gauss law at the boundary and the fact that net force 
experienced by an ion at the center of the sphere must be zero. 



Along with the above initial and boundary conditions, solution of SCF equations need 
to be obtained under additional constraints due to the fixed number of monomers (A^), ions 
{rij) and global electroneutrality so that 

ZpaN + Z^n^ = 0, Z+n+ = -Z_n_ = 0.6023 c,fi, (29) 

where Cg is salt concentration in moles per liter (molarity) and Q is in units of nm^. 



C. Numerical Technique 

We have solved SCF equations (Eqs. ( [7]- [H])) in real space using an explicit finite 
difference scheme for Poisson-Boltzmann-like equation (Eq. ([IS])) and the standard Crank- 
Nicolson^^ scheme for solving modified diffusion equation (Eq. ( [E])). Due to the use of 
properly normalized equations for densities, all the constraints mentioned earlier are always 
satisfied during the computation. 

As the solution of SCF equations is invariant when an arbitrary constant is added to the 
fields, so this constant needs to be fixed in order to obtain a unique solution. Choice of 
fixing this constant depends on the numerical scheme used in solving the SCF equations. 
We simply choose ip{R) = and J drri{r) = 0. We must point out that any method of fixing 
this constant does not affect the densities but the free energy gets changed by a constant. 

To realize the constraint f drri{r) = 0, ^ J drrj^r) is subtracted from the computed ri{r) 
at each iteration. This procedure leads to J drri{r) = in the final solution. Also, all the 
integrals are evaluated by the standard Trapezoidal^ rule and Broyden's method^^ is used 
to solve the set of non-linear equations. 



D. Reference System 



The choice of a reference system in free energy calculations depends on the physical quan- 
tity of interest. One of our goals in this study is to investigate the role of the polyelectrolyte 
chain in the free energy of the system. To study the role of the chain, the spherical cavity 
without any polymer (with small ions and solvent inside) is the appropriate choice. 

In the absence of the chain, free energy becomes 

F{pp = 0}-Fo{N = 0} ^ r. n, 
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(30) 



Note that in the absence of the chain, Uglt' = Q, due to the incompressibihty condition. 



Ahhough we have computed the free energy of a confined chain, it is worthwhile to 
compute a physical observable, which is readily measurable experimentally. So, we have 
computed the osmotic pressure of a confined chain by carrying out the variation of free 
energy (Eq. ( fTSl) ) with respect to the number of solvent molecules, but keeping the number 
of monomers and salt ions fixed ( and taking care of the fact that volume has to be changed 
to alter the number of solvent molecules for the incompressible system under investigation 
here). The osmotic pressure is given by^ 



where j = c, +, — . Within the saddle point approximation, F = F* and {F* — FQ)/kBT = 
—Hg for pure solvent (i.e. when N^rij = 0). Using the above formula, osmotic pressure 
comes out to be (within radial symmetry) 



where pe{R~) is the total charge density (in units of electronic charge) at a point close to the 
surface of the cavity and ^ in the superscript depicts the fact that saddle point approximation 
for the free energy has been used in deriving the result. Due to the coarse grained model 
used in studying the single chain, information at a length scale below Kuhn's segment length 
is not correctly captured by the model. So, R~ represents the point, which is at a distance of 
one Kuhn segment length from the surface. This point is well-discussed in the literature^"^ 
and will not be pursued further. 

If we were to imagine a system, where densities and fields are constant (as in the concen- 
trated bulk system), then the above expression simplifies to the well-known^^ expression for 
the osmotic pressure of a homogeneous system (represented by 11^) 



E. Osmotic Pressure: Contact Value Theorem 




(31) 




(32) 




(33) 



where Pp = Nb^/Q. In the hterature^^, osmotic pressure for fluids near surfaces is given by 
the fluid densities near surfaces (the so called "contact value theorem") and Eq. ([321) 
nothing but the analog of the "contact value theorem" for the inhomogeneous system with 
interactions. 



F. Mean Activity Coefficient 



Using the saddle point approximation for the free energy, we have computed the electro- 
static chemical potential of small ions. Carrying out the variation of free energy (Eq. ([IE])) 
with respect to number of small ions, the mean field estimate for the electrostatic chemical 
potentiali^ of small ions comes out to be 
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for j = c,+,-. (34) 



This is a straightforward generalization of the electrostatic chemical potential for homo- 
geneous system to an inhomogeneous one. Using these expressions, chemical potential for 
so that Vj^ = v_ = 1 for the monovalent salt) can be written as 



the salt Ay^By 



(35) 



However, individual activity coefficients or the chemical potential can not be measured 
experimentally. So, we construct the mean activity coefficient (7±) for the binary saltr^ 
defined by 
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where i/ = z/+ + z/_. Using Eq. ( 1351) and ( 1361) . mean activity coefficient is given by 



7± 



(37) 
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Using Schwarz's inequality^, it can be shown that the mean activity coefficient is always less 
than or equal to unity. For the spherical cavity with pure solvent and monovalent salt ions, 
potential is essentially constant everywhere (local electroneutrality) and that leads to mean 
activity coefficient being unity. Note that the Eq. ( 1371) can also be derived by considering 
a Donnan equilibrium between the interior containing polyelectrolyte chain and the exterior 
containing salt ions with salt concentration Cg. 



III. RESULTS 



Having presented the field theoretical treatment of a single flexible chain in the presence 
of solvent, we present the results obtained after solving SCFT equations for a negatively 
charged chain {Zp = —1) with monovalent salt. In this study, we have taken all small 
ions (counterions and co-ions) to be point charges and while calculating counterion density 
profiles, we have added contributions coming from the counterions released by the polyelec- 
trolyte chain and the salt. Also, it should be noticed that due to the point nature of these 
charges, counterions and co-ions are not excluded from the confining boundary, which can 
be done, in principle, by the introduction of an arbitrary wall potential^. Here, we simply 
assume that the depletion layer"^ for counterions and co-ions is very small as compared to 
the monomer due to an order of difference in their sizes and its negligible effects on the 
system properties. 

All the results reported in this paper have been obtained with a grid spacing of Ar = 0.1 
and chain contour discretization of As = 0.01 after putting h = Inm. 

1. Monomer and charge distribution 

To study the effect of confinement, we have varied degree of polymerization (A^) keep- 
ing all other parameters fixed. Increase in A^ with fixed R leads to a more confined 
environment with the increase in monomer density everywhere inside the cavity. This 
can be seen in Fig. [T|, where we have plotted monomer densities for a single poly- 
electrolyte chain for different values of A^ after choosing a particular set of parameters 
iXps = 0.45, a = 0.1, Is/b = 0.7, R/b = 5, = O.IM). These values for the parameters 
were chosen to mimic a salty flexible polyelectrolyte chain under spherical confinement in 
the presence of water as a solvent. Also, for comparison purposes, we have plotted monomer 
densities for the corresponding neutral chains in an athermal solvent {xps = 0) and in an 
equivalent solvent condition {xps = 0.45). Comparing the monomer density for the neutral 
chain in athermal solvent with a poorer solvent, it is clear that the exclusion of solvent 
molecules from the core of the coil is stronger as the solvent quality is decreased, as ex- 
pected. Now, when the polymer in the less good solvent is charged, the solvent exclusion 
effect is slightly weaker. This difference becomes negligible as the monomer volume frac- 



tion increases. Moreover, monomer density profiles for neutral and the polyelectrolyte chain 
in the same solvent are almost identical, which corroborates the coil conformation of the 
polyelectrolyte. 

We have also varied R by keeping the number of salt ions and a fixed, and the results 
obtained for radial densities are shown in Fig. [2] for a particular value of N (number of 
salt-ions corresponds to Cg = O.IM for R/b = 5). The observed difference in small ion 
density profiles for different -R's is attributed to the fact that total salt concentration (c^) 
is changing when R is being changed with the number of small ions kept fixed during this 
variation. In this figure, we have chosen = 100 and the rest of the parameters are the 
same as in Fig. [1] for the polyelectrolyte chain. From Fig. [2], it is evident that due to 
exclusion of the monomers from confining sphere surface, decrease in R leads to an increase 
in concentration of monomers in the interior. Had it been a bulk situation {R oo), small 
ion densities would have reached their bulk value (c^). Looking at small ions' densities in 
Fig. [21 it is clear that for the confined spherical system, this common presumed result is no 
more valid and small ion density profiles haven't reached their bulk value. In other words, 
net radial charge density is no more zero near the confining boundary (as shown in Fig. [3]) 
and a double layer systent^ is set up as a result of depletion of the chain from the surface. 

In this study, small ions are treated as point charges and hence, are not excluded from the 
surface of the spherical cavity in contrast to the monomer. As a result, when is increased 
keeping a fixed at stronger confinements (in terms of higher monomer densities), there are 
more number of counterions (positively charged) generated by the chain and hence, charge 
density near the surface of the cavity increases (Fig. [3]) . On the other hand, near the center 
of the cavity, the local charge density curves move toward zero as A^ is increased (i.e. interior 
of the cavity is becoming locally electro neutral). For lower degrees of confinements (i.e. low 
Pp), increase in A^ leads to increase in charge density everywhere due to more number of 
ions in the system. However, this behavior is seen over a very small density regime. 

We have also varied other parameters {cs,a and Ib)- Effects of all these parameters 
on monomer densities are consistent with the observation that a polyelectrolyte^ chain at 
high salt concentration is equivalent to a neutral chain with an effective excluded volume 
parameter given by 



1 -Pp 



1 




(1) 



So, the polyelectrolyte chain shows higher expansion^ as compared to its neutral analog^. 

On the other hand, at low salt concentrations, the chain can tend to attain a rod-like 
conformation due to dominance of electrostatic repulsions along the backbone and spherical 
symmetry is broken. In our mean field study, we use spherical symmetry and Gaussian 
model for the polyelectrolyte chain. Due to these limitations, spherical symmetry breaks 
down in the extremely low salt regime and we are not allowed to explore the low salt regime 
using the current model. 

These results show that confinement of the chain leads to the development of a charge 
density wave inside the cavity and an outcome of this charge density wave is that there is 
a potential difference across the center of the cavity and the surface. It is to be noted that 
for a confined system with only small ions inside, potential is constant everywhere and local 
electroneutrality is the equilibrium state (trivial solution of SCF equations). However, it is 
the depletion of the chain from the cavity surface, which leads to the accumulation of charges 
near the surface and as a result inhomogeneous charge distribution is attained. To study the 
role of polyelectrolyte in the free energy, we compare the free energies of the inhomogeneous 
phase with the cavity containing pure solvent and salt-ions in the next section. 



In Fig. m we have plotted different contributions to free energy for a salty system using 
Eqs. ([I9]- [23l) . Analyzing the contributions to free energy, it is clear that the free energy 
has four major contributions. At lower polymer volume fractions, free energy is dominated 
by entropy of small ions {—TSions/kBT) and solvent {— TS solvent /ksT). As volume fraction is 
increased, chain conformational entropy and polymer-solvent interaction energy also become 
important. Also, electrostatic energy part in free energy is small compared to other terms, 
stressing a minor role played by electrostatic energy in the crowded environment under 
investigation here. Share of each contribution to the free energy depends on the degree of 
confinement. For instance, when degree of confinement is extremely high {Nb'^/Q —>■ 1), 
solvent entropy and polymer- solvent interaction energy terms are negligible and free energy 



2. Free energy within saddle point approximation 



has two major contributions - small ion entropy and chain conformational entropy. Strictly 
speaking, our theoretical model breaks down as soon as monomer density inside the cavity 
is close to unity because the finite size of small ions and nature of interactions between the 
various species become important. All these effects can not be captured with our theory, 
which involves only two body interaction potential. 

A few comments about the shape of the plots in Fig. H] are due here. Shape of small ions 
entropy term is governed by number of small ions term (^^. J drpj{r) in Eq. ([22])) and when 

is increased keeping degree of ionization (a) and R fixed, number of counterions increases 
and hence, —TSions/ksT decreases with increase in A^ (almost linearly). Similarly, when A^ 
is increased while keeping R fixed, total number of solvent molecules in the cavity decreases 
(due to incompressibility) and hence, —TSsoivent/kBT increases. Increase in A^ for a fixed 
R/b leads to lower number of conformations available to the chain and hence, entropy of 
the chain decreases or —TSpoiy/ksT increases. Shape of excluded volume interaction energy 
and electrostatic interaction energy can be understood by the fact that in the asymptotic 
limit with respect to degree of confinement, the product of densities (monomer and solvent) 
as well as the electric potential are small. 

In implicit solvent computations^^, free energy goes linearly with N'^/Q. In order to see 
whether the same linear law is followed by free energy in explicit solvent model, we have 
plotted free energy as a function of N'^/Q for different values of R/b in Fig. It is found 
that free energy follows the linear law only for lower values of pp. For higher pp, there are 
deviations from this linear law due to the conformational entropy of the chain. Overall, the 
shape of the free energy curve is of the form Ylj=c + - s^j [I'^l^j/^) ^ -'-] lower pp and 
systematic deviations are seen for higher pp. 

Also, to highlight the role of the polyelectrolyte in the free energy of confinement, we 
have plotted the difference between the free energy of the spherical cavity with and without 
chain (i.e. AF* = F* — F{pp = 0}) for different radii of the confining cavity and monomer 
densities in Fig. [6l In these computations, it has been assumed that Wpp = Wss to get 
rid of an uninteresting constant. The results reveal that the confinement of the chain is 
a thermodynamically unfavorable process and for the same density, larger value of AF* 
for larger spherical cavity can be attributed to the difference in conformational entropy and 
polymer-solvent interaction energy for chains of different lengths in spheres of different radii. 

In scaling theories^, it is commonly asserted that the confinement free energy is extensive 



in A^. In order to confirm this assertion, we have plotted (F* — FQ)/NkBT for various values 
of monomer densities at different values of R (see Fig. [7|). Our results clearly show that the 
free energy under spherical confinement is not extensive in N, in contrast to the rectangular 
and cylindrical confinements^^. 

3. Osmotic pressure and mean activity coefficient 

In Fig. [HI we have plotted the osmotic pressure obtained from Eq. ( [32]) as a function 
of monomer density. For comparison purposes, we have also plotted the osmotic pressure of 
the homogeneous phase (Eq. ( [H^ ). It is found that the osmotic pressure is the same for 
the homogeneous and the inhomogeneous phase at lower monomer densities and ideal gas 
law for osmotic pressure is obtained in this regime. However, at higher densities, deviations 
from ideal gas law are seen and these deviations are larger for the inhomogeneous phase in 
comparison with the homogeneous phase. Comparing each term in Eqs. ( [321) and ( [33i) . it 
is found that the discrepancy between the inhomogeneous and homogeneous phase arises as 
a result of the depletion of the chain from the spherical surface. Due to the depletion, the 
monomer density at one Kuhn step away from the surface is higher for the inhomogeneous 
system in comparison with the homogenous system. Moreover, the log term involving solvent 
density and quadratic excluded volume interaction term involving the chi parameter add to 
the discrepancy. Other than these terms, electrostatic and small ions terms do not change 
much. Effect of the cavity radius on the osmotic pressure profiles can be easily explained 
using Eq. ( [HHl) and the fact that number densities of salt ions is higher for smaller cavity, 
when number of ions is kept fixed during the computation. 

In the absence of the polyelectrolyte, the local electroneutrality is the equilibrium state 
for the small ions. That means the mean activity coefficient for the salt is unity. However, 
due to the depletion of the chain from the surface, local electroneutrality gets broken and 
that leads to a deviation from unity. In Fig. [9l we have plotted the mean activity coefficient 
for the monovalent salt and it is clear that the presence of the polyelectrolyte chain leads to 
deviation from ideal behavior (local electroneutrality). The shape of the activity coefficient 
curves can be understood by the fact that local electroneutrality is attained for very large 
densities also, where number of small ions is large and almost uniformly distributed. Effect of 
R on these plots can be explained by the fact that the extent of inhomogeneity in electrostatic 



potential increases with an increase in R. 



4- One-loop fluctuation corrections : narrow depletion zone approximation 

The saddle point approximation used in computing free energy (Eq. ( [TSI) ) is valid 
when the number densities of small molecules and monomers is high (i.e. concentrated 
regime). In order to capture the role of fluctuations in the concentrated regime, we have 
expanded the integrand of functional integrals over fields around the saddle point solution 
up to quadratic terms in fields (Appendix-B) so that the functional integrals to be carried 
out are Gaussian (one-loop calculations). For low concentrations (dilute and semi-dilute 
regimes), this treatment breaks down and other techniques have to be employed. 

Even at the one-loop level, it is very difficult to sum the infinite series, which emerge as a 
result of Gaussian integrals . Moreover, these sums are plagued with ultraviolet divergences, 
which have to be regularized. However, we have been able to sum these series in the long 
chain limit (A^ — oo) when nR — ^ oo and S,^^R — > oo^ where and ^ are Debye and 
Edwards' screening lengths, respectively. In the concentrated regime for a very long chain, 
the width of the depletion zone near the surface of the sphere is very small as compared to 
the radius of the sphere and we can suppress the radial dependence of the densities (cf. Fig. 
[Hand E]). Taking this approximation and ignoring the correlation energy of the charges 
along the backbone of the chain (which is very small for a weakly charged polyelectrolyte), 
the free energy at one-loop level can be written as (Appendix-B) 



p* -j^ oo oo 



^ ^ k=l 1=0 



In 1 + + In 1 + 



(2) 



where = 4:7rlBj2c,+,-Zjnj/n,^-'^ = 12(y^ - 2xps)Nb/n and m is k^'' zero of the 
spherical Bessel function of order / (i.e. jiii^ki) = 0). The infinite sum over k can be 
computed exactly. However, the sum over / diverges. The divergence can be regularized by 
introducing an upper cutoff M on / and identifying cut-off independent part. The finite sum 
over / has been computed in Ref.— . The cut-off independent part gives 



+ ^''^^ + ^ (1 + 2 In «:/?) + ^ (1 + 2 In C'R) + ^R, 



ksT ksT 12n 327r ' ' 32n ^ ^ ^ 3 

(3) 



where S is the surface area of the sphere (= AttR'^). In the above expression, cubic terms 
in K and represent the bulk contribution of fluctuations to the free energy and other 
terms correspond to the finite size corrections, which arise as a result of confinement. Note 
that for R oo, these finite size contributions vanish and the well-known screening result 
is obtained. Also, the finite size contributions are smaller than the bulk contributions in the 
limit discussed here and hence, overall, fluctuations lower the total free energy. 

Fluctuation corrections to the osmotic pressure and the mean activity coefficient can be 
estimated in a straightforward way by using Eq. ( [S]). It must be kept in mind that Eq. 
( [3l) is strictly valid in the strong screening regime for an infinitely long chain. In order for 
these conditions to be realized, spherical cavity has to be very large to accommodate the 
long chain (because of incompressibility condition). On the other hand, numerical solution of 
SCF equations for large spherical cavities with a very long chain becomes very expensive. So, 
in this study, we limit ourselves to a qualitative discussion about the fluctuation corrections 
to osmotic pressure and mean activity coefficient. 

Qualitatively, the osmotic pressure is decreased due to fluctuations (because the fluctua- 
tions at one loop level lower the free energy). It can be shown quite easily that the leading 
corrections to the pressure profiles^*^ will be of the form —k,^/24tt and — ^^^/247r. In the 
thermodynamic limit of infinite volume for an infinitely long chain, these expressions become 
exact. As the osmotic pressure must be positive, this analysis sets the range of validity of 
the one-loop calculations^^ to be well above the overlap concentration (in the concentrated 
regime). In fact, below those concentrations, perturbative treatment to capture the role of 
fluctuations fails and non-perturbative methods to capture the role of higher order terms in 
the expansion of the integrand have to be used (e.g. in dilute and semi-dilute regime )^»^. 

The finite size corrections to these laws (originating from the last three terms on the right 
hand side in Eq. [H]) are small but positive and overall, increase the range of validity of the 
fiuctuation analysis. Similar analysis can be carried out for the mean-activity coefficient of 
the monovalent salt and the leading corrections coming from fluctuations^^ are of the form 
exp {—k,Ib/2). Of course, the analysis is valid for the low salt concentrations (Debye-Hiickel 
regime) because the size and the nature of the short-range interactions of salt ions become 
important at higher concentrations. 



IV. CONCLUSIONS 



We have studied a weakly charged flexible polyelectrolyte chain under spherical confine- 
ment using SCFT. In Sections [IT] and IIIIl we have demonstrated that SCFT predicts 
creation of a charge density wave and a potential difference across the center of the sphere 
and the boundary. 

We have also shown that for a given charge density along the backbone, free energy of 
a flexible chain {F* — Fq) has four major contributions - entropy of small ions, entropy of 
solvent , energy due to polymer- solvent interactions and conformational entropy of the chain. 
Share of each contribution to the free energy depends on the monomer density, degree of 
ionization and salt-concentration inside the sphere. However, electrostatic interaction energy 
plays a minor role in free energy for the weakly charged flexible polyelectrolyte. Our results 
show that the free energy is not extensive in the number of monomers. 

Osmotic pressure for the polyelectrolyte chain follows ideal gas law in the low monomer 
density regime and substantial deviations are seen as the monomer density is increased. 
Mean activity coefficient for the monovalent salt show a small, yet systematic deviation 
from unity, highlighting the role of the polyelectrolyte in breaking the local electroneutrality 
condition seen in the absence of the polyelectrolyte. 

One-loop fiuctuation analysis (within the approximation of narrow depletion zone) reveals 
that fluctuations lower the free energy and the free energy has additional contributions com- 
ing from screening effects due to monomers as well as small-ions. Without any confinement, 
the fiuctuations corrections at one-loop level^ have been shown to lower the free energy 
and the range of validity of the fiuctuation analysis is set by the condition that osmotic 
pressure must be positive. The concentrations over which the condition is satisfied, come 
out to be well-above the overlap concentration in polymer solution theories. As the finite 
size corrections to the fiuctuation contributions in the free energy are positive, so the range 
of validity of the fiuctuation analysis gets widened due to the confinement. 

Finally, we comment on the assumptions used in arriving at the above conclusions: 
(1) We have used the spherical symmetry of the system which is valid as long as coil 
conformation of the chain is retained and a rod like conformation is avoided. Also, we have 
used SCFT for a fiexible polyelectrolyte chain, which means the chain must not be stiffened 
due to charges on the backbone and large torsional barriers. These criteria are realized in 



the presence of high salt, where electrostatic interactions get screened and become short 
ranged. 

(2) Saddle-point approximation breaks down in the extremely dilute limit for monomers, 
solvent and small-ions. The approximation is strictly valid in the concentrated regime in 
the presence of enough salt and solvent. The necessity to have many solvent molecules, sets 
upper limit to be away from extremely dense regime, where Nb^/Q —>■ 1. 

In this study, fluctuations in densities and fields about the mean field solution of SCF 
equations have been treated perturbatively by expanding free energy functional up to 
quadratic order about the saddle point. This fluctuation analysis breaks down for low 
density regime, where higher order terms also contribute and need to be taken into account 
using an appropriate non-perturbative treatment. 

(3) Role of solvent is taken into account by taking the volume of a solvent molecule as the 
same as that of a monomer. In reality, solvent size may be smaller than Kuhn step length 
and incompressibility constraint is violated near the boundary. For the present study, we 
have simply assumed that the effect of the depletion zone for solvent on the system properties 
is negligible and can be ignored. All these limitations of the current model can be removed 
by taking the solvent size to be different from the monomer— and using Pp{r) + Ps(r) = p(r), 
where p(r) is a suitable function, which is po away from the boundary and falls from po to 
zero in a smooth fashion within the depletion zone^. However, the choice of p(r) and width 
of depletion zone is arbitrary and depends on the numerics of the problem. 

Similarly, small ions have been treated as point charges, which might not be a bad 
approximation knowing the fact that typical ion radii for monovalent cations and anions lie 
in the range ~ 0.06 — 0.22 nn>^. 

(4) In principle, the degree of ionization (a) should also be computed by minimization 
of free energy^ with respect to a. In order to simplify the numerical work, we have taken 
the degree of ionization (a) to be independent of Ib and avoided any ion condensation effects. 

(5) While splitting the free energy into energy and entropy, we have assumed that di- 
electric constant (e) of the solvent is insensitive to temperature. However, the temperature 



dependence of the dielectric constant^i^ can be easily incorporated in our theory. 
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APPENDIX A : SCFT with explicit solvent 

Here, we present a summary of the steps in obtaining the saddle point equations for a 
single polyelectrolyte chain in the presence of salt ions and solvent molecules (section (III Ap 



). The procedure is similar to the ones presented in Refs.— i^. As the first step, we define 
microscopic densities as 



1 



Nb 



p,(r) = - J ds6ir-Ris)) (A-1) 
Pji^) = ^^(^ - J = s,c, +,- (A-2) 

(A-3) 



i=l 



aZpPp{r) + ZjPj{r) 

i=c,+,- 



where pp{r), pj{r) and Pe(r) stand for monomer, small molecules (ions and solvent molecules) 
and local charge density, respectively. 

As the second step, we use a functional integral representation for the incompressibility 
constraint 

Y[5 (ppir) + psir) - po) = /"/}[^_^(r)]e-^/'^™+W('^''W+^=«-''«), (A-4) 

r 

where w+(r) is the well-known pressure field which enforces the incompressibility constraint 
at all points in the system and i = ^/—l. 

As the third step, dimensionless Flory's chi parameter is introduced to club together the 
three excluded volume parameters so that 

Xpsb' = wp, - !^!lE±^. (A-5) 



Note that the clubbing of excluded volume parameters has been possible because of the 
incompressibility constraint so that only one independent parameter appears in the theory. 
Otherwise, there would have been three independent parameters. 
Following these three steps, Eq. ( [T]) becomes 



exp 



-i j drw+ir) (/3p(r) + /3,(r) - po) - Xpsh^ J c?rpp(r)p,(r) 
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dr / dr 



, pe(r)pe(r') 
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eKBi |r — r " ' ' ^ ' 

So far, we have written the partition function in terms of the local densities. Now, we 
want to write this partition function in terms of the order parameter of the system. There 
are many different choices we can make for the order parameter— i^. However, in this study, 
we follow the method as described in ref.— to introduce the order parameter. As the fourth 
step, we define 



P+W = Pp(r) + p,(r) 
P-W = Pp(r) - ps{r) 



Po 



(A-7) 
(A-8) 



where p_(r) is the order parameter for the inhomogeneous system. 

As the fifth step, to go from densities to fields, we use the Hubbard- Stratonovich trans- 
formation for short range excluded volume interactions as well as long range electrostatic 
interactions, so that 
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Note that, w+(r), w_(r) and ip{r) are the real fields, which can be envisioned as chemical 
potential fields^!. Now, taking these two steps, Eq. ( IA-6P becomes 

exp(^--^^ = -j^j D[w+{r)] j D[w^{j)] j /^[^(r)] exp [-/{«;+, ^}] , 

(A-13) 



where 



f {w+,w-,il)} = f drw'i{r)-ipo ( (irw+(r) - InQp - ^ n^ lnQ 

j dr^(r) Vr ^(r) + y (^iV^i;^^ + + ^fipo 
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D[R(t)]exp -- / dtl-{ — \ +zw+{R}-«;_{R} + «Zpa^{R} 

(A-15) 

= j dr exp[-{iM;+(r) +w_(r)}] (A-16) 
= J dr exp[—iZj-ip{r)] for j = c, +,— . (A-17) 

The functional integrations over real fields are to be carried out by contour integration 
techniques and are almost impossible to compute exactly. However, if the number of small 
molecules is large then we can compute the integrals in the numerator by steepest descent 
technique, using the knowledge that the integral along the constant phase (imaginary part 
of the integrand) contour is dominated by the local minima of the integrand^ (note that the 
functional integrals in the denominator are divergent and ignored at the level of saddle point 
approximation). So, we approximate the functional integrals in the numerator by the value 
of the integrand at the local minima (where the phase comes out to be zero) so that the final 
approximated integral is real ( = / [w\,w'^,ip*} )• However, the saddle point values for w+ 
(= w'^) and ip{= come out to be purely imaginary in contrast to w_(= u^*), which is 
real. 

Also, it should be noted that the densities remain unchanged with the shift in fields by 
an arbitrary constant. So, we write iw-^-{r) — ^Xpsb^Po = iw^i^r) to get rid of the constant 



in saddle point equations. Now, deriving local minima equations for the integrand in Eq. 
( IA-14P with respect to w_ and ip, and using notation iw^{r) — w_(r) </'p(r), iw+{r) + 
w^{r) 0s (r), iip{r) — > ^/'(r), iw+{r) — > //(r), Eqs. ( [3- [I3]) are obtained. Using these 
saddle point equations and employing Stirling's approximation for Inn!, we obtain Eq. ( [TSll . 



APPENDIX B : Fluctuations around the saddle point 

Here, we provide details for one loop treatment of the fluctuations^. We expand the 
integrand / in Eq. ( IA-131) up to second degree terms in w+, W- and ip around their respective 
saddle point values. For convenience in writing the expansion, we introduce a dummy 
functional variable Cp{r), where p = 1, 2, 3 correspond to w+, w_ and ip, respectively. In this 
notation, the expression for free energy becomes (cf. Eq. ( lA-lSp ) 
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Kmp{r,v') 



^V{Cl,C2,C3} I 



(B-2) 



In the expansion, linear terms in fields vanish because of the saddle point condition. Kmp 
can be computed in a straightforward way and are presented here for completeness 
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In the case of confined chain within a spherical cavity, where electrostatic potential and 
monomer density are zero at the boundary (Dirichlet boundary conditions), the functional 
integral over (p can be carried out in Eq. ( IB- 11) by expanding all vectorial quantities in 
terms of spherical harmonics {Yim) and spherical Bessel functions {ji). For an arbitrary 
function h (r, r ') this means 

oo oo I 

h{r,r') = Y,Y.Y1 hMdki{r)dM{r')YiU0A)YLiO'A') (B-13) 
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where Uki is k^^ zero of the spherical Bessel function of order / (i.e. jiiiyki) = 0). Similarly, 
we expand (p{r) as 

oo oo I 

(pi^) = Y.Y.Y1 Cpkdki{r)YiUO,<P)- (B-15) 

k=l 1=0 m=-l 

Now, using orthogonal properties of spherical Bessel and spherical harmonics, the functional 
integrals can be computed. After some lengthy algebra, 
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In the above equation, the second term on the r.h.s. represents the contribution of 
correlations arising due to excluded volume interactions between monomers (neutral polymer 
contribution), the third term arises as a result of small-ions density fluctuations and the last 
term correctly represents the correlation energy of charges along the backbone of the chain. 



Also, C^/ and L^^^ correspond to the coefficients in the spherical harmonics expansion for 
the inverse operator of C (r, r ') and L (r, r '), respectively. Inverse operator for any arbitrary 
functional h (r, r ') is defined by 

j dv'h{r,r')h-\r\r") = 5(r-r"), (B-17) 

which gives Yl'^=i ^kuh~i = Sm, where 5ki is Kronecker delta. Also, L (r, r ') is given by 

L(r,r') = A(r,r') + 5(r,r') +4 J dr" J dv'" B {v,v") {C - AB)-\v\t"') B {v'\v') . 

(B-18) 

Unfortunately, computation of the coefficients involved in Eq. ( IB-161) requires three 
dimensional calculations for the densities and it is very hard to compute the sums exactly, 
which diverge in general. To gain an insight into the problem, we have identified one 
particular case, where we can evaluate the first two terms analytically. 

If the number of small molecules (solvent and small ions) is large, then operators Bj (r, r ') 
become diagonal and also, number densities (pj(r)) show a weak dependence on r. Similarly, 
if degree of polymerization is large (strictly if A^ — oo), then monomer density becomes 
independent of r except near the surface of the cavity characterized by the width of the 
depletion zone, which can be neglected. If both of these conditions are satisfied, then 
suppressing position dependence of densities, we get 

^oooo / oooo X^oooo 
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where pp = and Wr = '^Xps is the renormalized excluded volume interaction 

parameter. For constant monomer density (while working in the limit A^ — > oo), A{r,r') 
satisfies 

Vr'^(r,rO = -^5(r-r') (B-20) 

so that b'^Aki = 12NhB? / {vIiVL) . Neglecting the correlation energy of the charges along the 
backbone of the chain (the last term on r.h.s. of Eq. ( IB-16P ) and an ultraviolet divergent 
sum in Eq. ( IB- 191) (ffist term on r.h.s. of the equation), Eq. ([2]) is obtained. 
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FIGURE CAPTION 

Fig. 1.: Effect of on monomer densities - comparison witli the corresponding neutral 
chains . In above plots, we have chosen R/h = 5 and N = 100, N = 200 and N = 300 
from bottom to top, respectively . 

Fig. 2.: Effect of cavity radius (R) on ion densities. In these plots, we have chosen = 
100, a = 0.1, Is/b = 0.7, Xps = 0.45 and number of salt ions is kept fixed (in all 
these plots, number of salt ions is equivalent to salt concentration of O.IM for a 
sphere of radius R/b = 5). Solid, dashed and dash-dotted lines represent p(r) = 
app{r), Pcij) + p+(r) and P-{r), respectively. 

Fig. 3.: Effect of cavity radius (i?) on net charge density, pe{r) = J2j=c + - ^jPji^) + 
Zp0tpp{r). In these plots, we have chosen a = 0.1, Is/b = 0.7, Xps = 0.45 and number 
of salt ions is kept fixed (in all these plots, number of salt ions is equivalent to salt 
cone of O.IM for sphere of radius R/b = 5). is increased in steps of 50 starting from 
50. 

Fig. 4.: Different contributions to the free energy within the saddle point approximation 
for salty systems. Here, we have chosen = 100, a = 0.1, Is/b = 0.7, Xps = 0.45, Cs = 
O.IM, R/b = 5. 

Fig. 5.: Effect of confinement on the free energy within the saddle point approximation. 
Here, we have chosen a = 0.1, Is/b = 0.7, Xps = 0.45 and number of salt ions is kept 
fixed so that the number of salt ions is equivalent to salt concentration of O.IM for 
sphere of radius R/b = 5. 

Fig. 6.: Difference in free energy of the spherical cavity with and without polyelectrolyte 
chain, AF* = F* — F{pp = 0}. All other parameters are the same as in Fig. [51 

Fig. 7.: Non-extensive nature of free energy of confinement. Parameters are the same as 
in Fig. [51 

Fig. 8.: Comparison between osmotic pressure for the confined polyelectrolyte chain and 
the homogeneous phase (cf. Eq. ([331)). Solid lines correspond to the inhomogeneous 
case and dashed lines represent the homogeneous system. Parameters are the same as 
in Fig. [SI 



9.: Mean activity coefficients for monovalent salt as a function of monomer density of 
the polyelectrolyte. For comparison purposes, parameters have been chosen to be the 
same as in Fig. [5l 
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